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the process. We show that, using a careful definition of equilibrium distribution and 
taken into account the appropriate time reversal transformation, usual equilibrium 
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of a covariant supersymmetric formulation of a multiplicative Markovian white-noise 
process and study some of the constraints it imposes on correlation functions using 
Ward-Takahashi identities. 
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1. Introduction 

From the early studies of Brownian motion, more than a hundred years ago p], the use 
of stochastic differential equations to model a wide variety of dynamical systems has 
grown dramatically. Applications can be found along a wide area of scientific research, 
from physics and chemistry [21 [3] , through biology and ecology (U [5] , to economy and 
social sciences [HI [7]. 

In the original formulation of Brownian motion, the influence of the medium on 
a diffusive particle is modeled by splitting its effect in two parts: a deterministic one, 
given by an homogeneous viscous force, and a stochastic part, given by a random force 
with zero expectation value. In this way, fluctuations exhibit as an additive noise and, 
consequently, the considered model represents an additive stochastic process described 
by a Langevin equation. However, the viscous force could have non-homogeneous 
contributions, for instance, in the presence of boundary conditions, such as a diffusion of 
a Brownian particle near a wall [HI E] - If the diffusion is not homogeneous, fluctuations 
could depend on the state of the system and they can be regarded as the product of 
a random force and a function of the state variable. In that situation, multiplicative 
noise is defined and the model is known as a multiplicative stochastic process. Another 
interesting example of multiplicative noise is the stochastic Landau-Lifshitz-Gilbert 
equation [10], used to describe dynamics of classical magnetic moments of individual 
magnetic nano-particles. In this case, the noisy fluctuations of the magnetic field 
couple the magnetic moment in a multiplicative way. The classification of additive 
or multiplicative noise should be considered in the more general context of external and 
internal noise [21 [TTJ. In the former case, dissipation and noise can be considered as 
two effects with different microscopic origin. In other words, in the absence of noise, 
the classical deterministic system is perfectly well defined. However, in the latter case, 
dissipation and fluctuation have the same intrinsic origin and it is not possible to "turn 
off" one of these effects. 

The theory of stochastic evolution provides a beautiful connection between 
dynamics and statistical physics. In general, Einstein relation or, more generally, 
fluctuation-dissipation relations associate dynamical properties of the system with 
thermodynamical equilibrium. However, stochastic dynamics not necessarily model 
physical systems, where the long time evolution should conduce to thermodynamical 
equilibrium. In fact, it is possible to have more general stationary state distributions 
that represent an equilibrium state in the stochastic dynamical sense, not related 
to thermodynamics. Moreover, stochastic dynamics provide an interesting approach 
to out-of-equilibrium statistical mechanics [121 HB]- Interestingly, it is possible to 
attribute thermodynamical concepts like heat, entropy or free energy to each trajectory 
of an stochastic evolution, given rise to the research field usually called stochastic 
thermodynamics [T3] , 

In this work, we would like to present a study of equilibrium properties of Markovian 
multiplicative white-noise processes. A prototype of these processes is represented by 
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the Langevin equation ([T]). It is well known that, due to the Gaussian white-noise 
distribution, the continuum limit of the discretized time evolution of multiplicative 
processes is not unique. In fact, there are different prescriptions to perform this limit. 
Perhaps, the most popular prescriptions are the Ito [15] and the Stratonovich [16] ones, 
each one producing a different stochastic evolution and forcing different rules of calculus. 
We will discuss in detail equilibrium properties in the more general prescription called 
Generalized Stratonovich convention [17] (also called "a-convention" [18] in the field 
theory literature). In this prescription, a continuum parameter < a < 1 is defined in 
such a way that each value of the parameter corresponds to a different discretization 
rule of the stochastic differential equation. The value a = corresponds with the Ito 
prescription while a = 1/2 is the Stratonovich one. 

The concept of equilibrium is tightly related with the concept of time reversal 
symmetry. For this reason, it is essential to carefully define the time reversal 
transformation of the stochastic process. In multiplicative white-noise processes, the 
forward and backward stochastic trajectories evolve with different prescriptions (except 
when considering Stratonovich convention). For instance, if Ito convention (a = 0) 
is defined for the forward evolution, the backward trajectory evolves with the so 
called Hanggi-Klimontovich [T7J [191 1201 [21] prescription (a = 1). In general, we will 
show that the prescription a is the time reversal conjugate of (1 — a). In addition, 
since the asymptotic stationary probability distribution also depends on the chosen 
prescription, the correct definition of time reversal transformation, compatible with a 
unique equilibrium distribution, is quite involved. 

It is very useful to use, instead of the Langevin or Fokker-Planck approach, a path 
integral formalism [22J, in which the central object of the theory is a set of stochastic 
trajectories. Using this formalism, we will analyze equilibrium properties, such as 
detailed balance relations, microscopic reversibility and entropy production. We will 
show that, due to a careful definition of equilibrium distribution and taken into account 
the appropriate time reversal transformation, usual equilibrium properties are satisfied 
for any value of a. Even thought the path integral representation of the stochastic 
process is suitable for the computation of correlation functions and responses, specific 
calculations are very cumbersome, since each value of a defines different differentiation 
and integration rules. In particular, the "chain rule" to compute total derivatives of a 
stochastic process depends on a. We make a generalization of the Ito formula for any 
value of the discretization parameter. In this work, we explicitly show the importance 
of this "generalized chain rule" , in the path integral formalism. 

The stochastic process can also be represented as a path integral in an extended 
functional space, by introducing auxiliary commutative as well as anti-commutative 
Grassman variables. The main advantage of this extension for multiplicative noisy 
systems resides in the prescription independent character of the formulation. While 
in the original path integral formulation the prescription appears as a continuum 
limit ambiguity, in the functional Grassman formalism, the ambiguity appears in the 
definition of equal-time Grassman Green functions [22] . In this way, provided we do not 
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integrate the Grassman variables, we can perform any calculation without specifying a 
particular prescription. In some sense, all the difficulty introduced by the generalized 
Ito calculus is being codified in a simpler Grassman algebra. 

As a by-product, it is possible to study hidden symmetries of the stochastic 
process. Indeed, the class of systems we have studied in this paper is invariant 
under linear transformations in the extended functional space. This hidden symmetry, 
called supersymmetry (SUSY), was recognized in other stochastic processes several 
years ago [231 EH EE]- SUSY properties have been extensively studied for additive 
stochastic processes [261 EI] arid for non-Markovian multiplicative processes |28j. In 
another hand, for Markovian multiplicative white-noise systems, important progresses 
have recently been reported [29]. The main difficulty in the SUSY formulation of a 
Markov multiplicative process is related with the great variety of prescriptions available 
to define the Wiener integral, which produces several stochastic evolutions with different 
final steady states. Moreover, time reversal transformations mix different prescriptions. 

In this paper, we deduce a covariant supersymmetric formulation of a multiplicative 
Markov process showing in detail the importance of the prescription dependent 
equilibrium distribution and the correct definition of time reversal transformation. 

Physically, SUSY encodes equilibrium properties of the system. Some of the 
constraints it imposes on correlation functions (Ward-Takahashi identities) are related to 
fluctuation-dissipation theorems [30]. This property has acquired a refreshing interest 
due to the growing importance of stochastic out-of-equilibrium systems [31]. In this 
sense, it is possible to understand out-of-equilibrium dynamics as a symmetry breaking 
mechanism. We analyze in detail Ward-Takahashi identities specially related with 
fluctuation-dissipation theorems. We show that this theorem is fulfilled independently of 
the chosen prescription to define the process. However, the linear response, fluctuations, 
as well as the equilibrium distribution, do depend on the specific prescription. 

The paper is organized as follows: In Sect. [2] we define our model by means 
of a Langevin equation with multiplicative noise. In Sect. [3] we carefully study the 
problem of time reversal symmetry and the equilibrium probability distribution. With 
these results, we analyze properties such as detailed balance and entropy production in 
Sect. HI In Sect. El we present a covariant supersymmetric formulation of the stochastic 
process and we deduce the fluctuation-dissipation theorem for arbitrary prescriptions. 



Finally, we discuss our conclusions in Sect. El leading for Appendix A some details of 
the calculations. 

2. Multiplicative white— noise stochastic evolution 

The main purpose of this section is to define the model and to establish the concepts 
and the notation we use in the rest of the paper. For simplicity, we consider a single 
random variable x(t) satisfying a first order differential equation given by 

^ = f(x(t)) + g(x(t))((t), (1) 
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where ((t) is a Gaussian white noise, 

(C(t)) = , (C(t)C(t')) = 8(t-t'). (2) 

The drift term f(x) and the square root of the diffusion function g(x) are, in principle, 
arbitrary smooth functions of x(t). The only restrictive condition is that g(x) should 
be "invertible" , i. e., g(x) 7^ 0,Vx. 

As it is very well known, to completely define equation (JTJ), it is necessary to give 
sense to the ill-defined product g(x(t))((t), since ((t) is delta correlated. The problem 
can be easily understood looking at the integral 

J g(x(t)) C(t)dt = J g(x(t)) dW{t) , (3) 

where we have defined the Wiener process W(t) as ((t) = dW(t)/dt. By definition, the 
Riemann-Stieltjes integral is 

/ g(x(t)) dW(t) = kfj^-py - W{tj)) (4) 

J 3=1 

where Tj is taken in the interval [tj, tj + i] and the limit is taken in the sense of mean-square 
limit [3j. For a smooth measure W(t), the limit converges to a unique value, regardless 
the value of Tj. However, W(t) is not smooth, in fact, it is nowhere integrable. In 
any interval, white noise fluctuates an infinite number of times with infinite variance. 
Therefore, the value of the integral depends on the prescription for the choice of Tj. 
There are several prescriptions to define this integral that can be summarized in the so 
called "generalized Stratonovich prescription" jTTj or "a-prescription" [IB] , for which 
we choose 

g(x(Tj)) = g((l — a)x(tj) + ax(tj + \)) with < a < 1. (5) 

In this way, a = corresponds with the pre-point Ito interpretation and a = 1/2 
coincides with the (mid-point) Stratonovich one. Moreover, the post-point prescription, 
a = 1, is also known as the kinetic or the Hanggi-Klimontovich interpretation [I7J [19j 

US [21]. 

In principle, each particular choice of a fixes a different stochastic evolution. 
In many physical applications, a weakly colored Gaussian-Markov noise with a finite 
variance [32] is considered. In this case, there is no problem with the interpretation of 
equation (CD) and we can take the limit of infinite variance at the end of the calculations. 
This regularization procedure is equivalent to the Stratonovich interpretation, a = 
1/2 [2j [26]. However, in other applications, like chemical Langevin equations [2] or 
econometric problems [6l [7], the noise can be considered principally white, since it could 
be a reduction of jump-like or Poisson like processes. In such cases, the Ito interpretation 
(a = 0) should be more suitable. Hence, the interpretation of equation (PQ) depends 
on the physics behind a particular application. Once the interpretation is fixed, the 
stochastic dynamics is unambiguously defined. Summarizing, in order to completely 
define the stochastic process described by equation ([I]) we need to fix a couple of 
functions (f,g) and the parameter < a < 1. 
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3. Equilibrium distribution and time reversal 

The a-prescription introduces some difficulties in the proper definition of time reversal 
evolution that we describe in this section. 

3.1. Fokker-Planck approach 

Using equation ([[]) and the a-prescription, it is immediate to obtain [33] a Fokker-Planck 
equation for the probability distribution P(x,t): 

where g' = dg/dx. This equation can be cast into a continuity equation, 

dP{x,t) + gjfot) = Q 
dt dx 
whith the probability current given by 

J(x,t) = [f(x) - (1 - a)g(x)g'(x)) P(x,t) - l -g\x)^^. (8) 

Then, the probability distribution P(x,t) is given by the solution of the Fokker-Planck 
equation (jHJ) supplemented with the initial condition P(x, 0) = Pm(x) and two boundary 
conditions. We can fix, for instance, the values of J(x, t) at two boundary points (xi,Xf), 
where eventually they can be taken as ±oo. However, due to the fact that P(x, t) is 
normalized to one, for any value of t, then J(xi,t) = J(xf,t). 

We suppose that, at long times, the probability rapidly converges to a steady state 
P s (x), given by 

P s (x) = lim P(x,t) = N e~ u{x \ (9) 

with the normalization constant iV -1 = J^dx e~ u ( x >. In this state, the stationary 
current gets the form 

J S {x) = Ne~ u ^ - (1 - a)g(x)g'(x) + ^(x)^^ (10) 

and the stationary Fokker-Planck equation acquires the simpler form 

^ = 0, (11) 

dx 

or, simply, J s (x) = J = constant. 

Thus, for a given value of J, we could have a stationary state characterized by 
the function Uj(x). The stationary probability density Pj(x) = Nj e~ Uj ^ does not 
represent, in general, an equilibrium state. An equilibrium state is characterized by 
a net zero stationary current (J = 0). When J / 0, although the current is already 
conserved, there is a stationary probability flux, characterizing an out-of-equilibrium 
regime. In most physical applications, the diffusion function g(x) is a polynomially 
growing function for large x. In these cases (and for a single variable), the only possible 
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solution is the equilibrium one P eq = Pq(x) = N e~ Uo<yX \ However, for exponentially 
growing functions, out-of-equilibrium steady state solutions are possible. 

Let us focus on the equilibrium state, defined as the solution of the stationary 
Fokker-Planck equation with zero current probability. Thus, for J = 0, there is an 
obvious solution of equations (ITDjl and (TTTj) . given by 

U cq (x) = -2 f 4rl dx ' + (!"«) ^g\x). (12) 
J g [x ) 

Therefore, given the functions f(x) and g(x), supplemented by a value of < a < 1, 
the probability density tends at long times to the equilibrium distribution 

P cq (x) = Ne~ u ^ x \ (13) 

where U eq (x) is given by equation f[T2"l) . provided e~ Uc ^ is integrable to define the 
normalization constant N. 

In the case that equation ([1]) models a conservative physical system, i. e. , if the 
otherwise deterministic system is characterized by an energy potential V(x), the drift 
function f(x) does not depend explicitly on time and is given by 

m = -\^) dJ ^- (14) 

This equation relates fluctuations (given by g(x)) with dissipation, therefore it can be 
considered as a local generalization of Einstein relation for Brownian motion. Using this 
expression, the equilibrium "potential" has the form: 

U cq (x) = V(x) + (l-a)\ng 2 (x). (15) 

As expected, the equilibrium distribution depends, not only on the given functions 
(f(x),g(x)) or, alternatively, on (V(x), g(x)), but also on the value of the a-prescription 
which defines the Wiener integral. If equation ([ID is used to model a physical 
system, which is expected to asymptotically converge to thermodynamic equilibrium 
(U eq (x) = V(x)), the only possible choice is the Hanggi-Klimontovich interpretation 
a = 1; for this reason, the election of a = 1 is sometimes called "thermal prescription". 
Recently, an experimental evaluation of a was reported [3U [35] in the study of a falling 
colloidal particle near a wall. In this work, the value of a = 1 was recognized as the 
proper one to make the correctly identification of the microscopic forces acting on the 
colloidal particle. With any other choice, the equilibrium potential U eq (x) ^ V(x), 
including the usual Ito (a = 0) and Stratonovich (a = 1/2) prescriptions. 

Except for the Stratonovich prescription a = 1/2, each value of a is associated 
with different calculus rules. For this reason, in most physical literature, Stratonovich 
prescription is preferred, because it enables "normal" rules of derivation and integration. 
Despite the Hanggi-Klimontovich convention being the only one which guarantees 
the convergence to the termodynamic equilibrium, calculations can be particularly 
cumbersome in this interpretation. However, if we prefer to work in any other different 
fixed prescription, and we insist in modeling a physical system which converges to 
Boltzmann's equilibrium, we need to modify equation ([T|). It is done by adding a 
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"spurious drift force" [33] which cancels the last term of equation ( TT5]) . That is, if we 
define a new a-dependent "potential" V a = V — (1 — a) lng 2 (x), then, the equilibrium 
distribution U eq (x) = V(x), for any value of a. 

In this work, we will not adopt this point of view. We will consider a stochastic 
process completely defined by fixing (f(x),g(x),a) in eq. ([I]). Thus, the probability 
distribution P(x, t) flows, at long times, to an a-dependent equilibrium distribution 
given by equation ( jl2j) or alternatively equation f|T5|) . In this way, we are able to study, 
in a unique formalism, model systems with general equilibrium distributions that go 
from Boltzmann thermal equilibrium to power-law distributions. A simple example of 
the latter case is to consider a pure noisy system with V(x) = 0. In that case, the 
equilibrium distribution is P eq ~ 1/ g 2 ^ l ~ a \x). 

We will show that this set of distributions, characterized by asymptotically zero 
current probability, is consistent with all the usual properties of equilibrium evolution, 
like, for instance, detailed balance and entropy production. 



3.2. Time reversal 

Consider, for instance, the Langevin equation ([I]), in the simplest case of f(x) = 0. 
Naively, we could compute the backward evolution of the stochastic variable x(t) by 
just changing the sign on the velocity x — > —x, in such a way that, for a given noise 
configuration, the particle turns back on its own feet. However, the situation is not so 
simple. Let us look more closely. Consider a time interval (t,t + At). The forward 
evolution is obtained by integrating equation ([I]) between the initial and final times t 
and t + At respectively, 

/•t+At rl r (t') rt+At 

x(t + At) - x(t) = jf °^-df = jf g(x(t>))at>)dt>, (16) 

while the backward evolution is simply obtained by changing the initial and final 
integration limits, 

x(t) - x(t + At) = f ^p-dt' = f g(x(t'))at')dt', (17) 
Jt+At at Jt+At 

where we use the notation x, just to differentiate backward from forward evolution. If 

the integrals were "normal" integrals, we could use the trivial property J a b = — J fe a . In 

that case, there would be no need to differentiate backward and forward variables, since 

equations f fT6|) and ffT7|) would be the same equation and x = x. However, the integrals 

are Wiener integrals and need to be carefully defined. As we saw in the last section, 

/•t+At n 

/ g(x(t)) dW(t) = kjj^^y - W(tj)), (18) 

J t ■ 1 



with 



g(x(r j )) = g((l-a)x(t j ) + ax(t j+1 )), < a < 1. (19) 
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The time reversal integral is obtained by changing tj in equation (ITS]) . The 

important point is that g(x(rj)) also depends on (tj,tj + x). Then, 
r t n 

/ <?(*(*)) dUT(t) = " (2°) 

Jt+At 3=\ 

with 

j(a;(r 3 )) = ? (ra(t 3 ) + (l-aMy), < a < 1. (21) 

where was obtained from equation flT9|) . by replacing tj -H- tj + i or, equivalently, 
a — >■ (1 — a). Therefore, the time reversed stochastic evolution is characterized by the 
transformations 

x(t) -> x(—t) and a -<->• (1 - a). (22) 

In this sense, we say that the prescription (1 — a) is the time reversal conjugate 
(TRC) of a. Then, the post-point Hanggi-Klimontovich interpretation is the TRC of 
the pre-point Ito one, and vice versa. The only time reversal invariant prescription 
is the Stratonovich one, a = 1/2. This means that, except for the Stratonovich case, 
the backward and forward stochastic paths do not have the same end points. This is 
illustrated in figure ([TJ, where we compute a "time- loop" evolution. Consider we want 
to compute the evolution of the system starting at x(t), going forward a time interval 
At and, then, turning back in time the same interval —At. 

The forward path in figure ([TJ) can be computed from equation (1T6|) considering 
x(t) as an initial condition, 

ft+At 

x{t + At) = x(t) + J g(x(t'))({t')dt' . (23) 

Then, the backward path is computed from equation f[T7|) considering x(t + At) as an 
initial condition 

x(t) = x{t + At) + f g(x(t'))({t')dt'. (24) 

Jt+At 

Replacing equation ( [23|) into ( [24|) . 

rt+At ft 

A a x(t) = x(t) - x(t) = / g(x(t'))((t')dt' + / g(x(t'))((t')dt', (25) 

Jt Jt+At 

and using equations ( fT8|) and ([20|) we find 

n 

A a x(t) = ^J2\g(x{Tj)) - g(x(r,))] (W(t j+1 ) - Wfo)) . (26) 

3=1 

We see that, in general, A a x(t) ^ since the backward and forward evolutions develop 
with different dual prescriptions. In the Stratonovich case, a = 1/2, g(x) = g(x) and 
Ax 1/2 (t) = 0. 

This fact is reflected in the a-dependence of the equilibrium stationary state U eq (x) 
(equation ( TT2"]) or equation (EES}). While the stationary state in the forward evolution 
is Uf w (x) = V(x) + (1 — a) lng 2 (x), the stationary state reached by the backward one 
is given by Ub w (x) = V(x) + a\iag 2 (x), and (except for the Stratonovich prescription) 
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forward 




x(t + At) 



backward 



Figure 1. Sketch of a "time- loop" evolution of Langevin equation (TjQ). The forward 
evolution t — > t + At is performed with a prescription a, while the backward trajectory 
t + At — > t evolves with the time reversal conjugate prescription (1 — a). Except for 
the Stratonovich convention a = 1/2, Ax(t) a ^ 0. 



Uf w 7^ Ubw Thus, this is not the backward evolution we are interested in. The key 
point is to determine the backward stochastic process which, at long times, converges 
to the same equilibrium distribution as the forward one. To obtain it we write the time 
reversed Fokker-Planck equation as 

dP(x,t) 

dx 



1 



a 



)gg>}P( X ,t)-l d ~ 



(g 2 P(x,t) 



(27) 



dt dx ^ > 2 dx 2 

where we made the substitutions x(t) — >■ x(—t), a — > (1 — a) and / — > f in equation (JB]). 
In equation (|27|) . P(x,t) represents the backward probability distribution of the event 
(x, t). Now, we need to specify the drift f(x) which produces a backward evolution 
that asymptotically converges to the equilibrium distribution P eq = Ne~ Ueq , with 

, we solve 



Ne~ 



U eq (x) given by equation (fT5j) . Setting dP(x,t)/dt = and P eq 
equation (1271) for /, obtaining 

f(x) = f(x)-(l-2a) g{x)g'{x). (28) 

Therefore, the time reversal transformation T which makes physical sense, meaning 
that it produces a backward evolution which converges to a unique equilibrium 
distribution, is given by 

x(t) — > x(—t) 



r 



a 



f 



1 



a 



(29) 



/ - (1 - 2a) gg> 



Notice that, since 1— 2a is odd under the transformation a — > I— a, the time reversal 
operator defined in this way satisfies T 2 = /, as it should be. As we have previously 
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stressed, in the Stratonovich prescription f = f and a = 1 — a, so, the time reversal 
operator T simply corresponds to change x(t) — > x(—t). In any other prescription the 
definition of T is more involved, given by equation ( 1291) . Although we have defined the 
time reversal transformation to match the correct equilibrium distribution in the long 
time limit ±00, this transformation also gives the correct answer in the stationary out-of- 
equilibrium case J 7^ 0. In this case, the stationary forward and backward probabilities 
are different, and they are related by J — > — J. That is, Pj w (x) = P[^(x) as can be 
easily seen by applying the transformation ( 1291) to the Fokker-Planck equation ([7]) and 



4. Detailed balance, microscopic reversibility and entropy production 

In this section we check that the concept of a-dependent equilibrium and the time 
reversal transformation introduced in the previous sections are consistent with the usual 
concepts of detailed balance, microscopic reversibility and entropy production. To do 
this, we use the path integral formalism, which turns out to be very convenient for 
handling stochastic trajectories. 



4-1. Path Integral formulation 



The transition probability P(xf,tf\xi,ti) represents the conditional probability of 
finding the system in the state Xf at time tf provided the system was in the state 
Xi at time £j. In the path integral approach, it can be written as [22] 



P(xf, tf\xi, U) = J Vx det 1 (g) e~ 



■S[x] 



(30) 



with the boundary conditions x(ti 



x, 



and x(tf) = Xf and the "action" S[x] given by 



S 



dt 



2g* 



dx , 

Tt~ f + a " 



1 2 



(31) 



The probability for the time reversal process P(xi,tf\xf, £j) is obtained by applying 
the time reversal transformation T of equation ( |29i) to the forward probability, 
equations (|30|) and ( ]3Tj) . 



P(Xi , tf J Xf , tj 



Vx det~ l (g) e 



Six] 



(32) 



with the boundary conditions x(ti) — Xf and x(tf) = Xi and the time reversed "action" 
S[x] = TS[x]. 

In order to compute S, it is interesting to write S[x] in another form to make 
explicit its properties under the time reversal transformation. Expanding the integrand 
of equation (!3l!) and considering a conservative system / = — (g 2 /2)V (equation (fi~4l) ) 
we find 

2 , .2 



S 



dt 



1 j dx 



+ 



1 n 
2? 



-g 2 V + agg' 



a 
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ld[V + alng 2 } 
2 dx 



(33) 



The last term can be transformed in a total time derivative. However, for arbitrary 
values of a, we need to apply the proper chain rule, since usual calculations rules are 
valid only for the Stratonovich convention a = 1/2. In the general case, for an arbitrary 
differentiable function of a random variable Y(x(t)), the chain rule reads 



dY{x{t)) 
~dt 



dY dx 
dx dt 



(1 



2a) d 2 Y 2 
9 ■ 



2 a* - ■ (34) 

Clearly, for a = 1/2, equation ( I34p is the usual chain rule. For a = 0, this formula is 
known as the ltd formula. For the Hanggi-Klimontovich prescription, a = 1, the latter 
differentiation rule differs from the Ito formula just in the sign of the last term. This 
is because both prescriptions are time reversal conjugate. In fact, the second derivative 
term is odd under a — )■ (1 — a), as it should be, since all terms in (IMj) should be odd 
under time reversal. 

Applying equation ( 1341) to the last term of equation ( |33|) . we find 



where 



S 



S 



a \ng' 



+ 3, 



(35) 



dt I— — 



t\2 




-agg'V - -g 2 V"- 



r\2 



(36) 



The most important property which results from equations (I35p and (I3^|) is that S 
is time reversal invariant, S = TS. Therefore, by using equation fl35l) . the time reversed 
action S can be easily determined: 

*/ 



S = TS 



■Itv 



(2 -3a) \ng 2 ) ' +S 



(37) 



Thus, for an arbitrary a, the variation of the action S[x] under a time reversal 
transformation is just a total derivative term. Using equations ( 1351) and ( 1371) . it is 
immediate to show that 



with U eq (x] 



S - S = U cq (x f ) 
given by equation ( TT5l) . 



U e q(Xi 



(3* 



4-2. Detailed balance 

One of the consequences of the time reversal operator ( 1291) is that the conditional 
probabilities P ^ P. However, using equations (!30|) . ( 132|) . (|35|) and (!37|) . it is simple to 
show that 

P(xf,tf\xi,tj) = e -[y+(i- Q )in 9 2 ]^ - 3g . 

P(Xi,tf\xf,ti) 
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satisfying the detailed balance relation 

P cq (x,i) P(Xf,tf\Xi,ti) = P(Xf,ti\Xi,tf) PeqiXf), (40) 

with P eq {x) given by equations ( fTBl) and ( |T5i) . 

So, we have shown that forward and backward conditional probabilities, with the 
definition of equilibrium state given by the effective potential of equation (fl5|) . satisfy 
the detailed balance relation, equation ( 1401) . for arbitrary values of < a < 1. In the 
special case of a — 1/2, P — P, and equation fj40l is the usual detailed balance relation. 



4-3. Microscopic reversibility and entropy production 

It is not difficult to analyze the previous results from the perspective of stochastic ther- 
modynamics [36]. The detailed balance relation (14*0!) represents a connection between 
a transition probability and its time reversed one. Moreover, the transition probability 
was obtained by integrating over all possible stochastic trajectories. However, there 
is also a relation between the probabilities associated with each individual stochastic 
trajectory, that is usually called microscopic reversibility [12]. To show up this relation 
explicitly, we can assign to each stochastic trajectory x(t), with end points (xj,tj) and 
(xj,tf), a weight 

V{x[t)\xit u xft f ) = det'\g) e~ s ^ , (41) 

The heat dissipated into the environment Q, and thus, the increase of entropy in the 
medium As m associated with one specific trajectory, is given by [37] 

V(x(t)\xiti,x f tf) 



As m = PQ[x(t)] = In 



(42) 



V(x(t)\xitf,x f ti)_ 

Hence, for each stochastic trajectory beginning in a state with initial distribution p(xi) 
and ending with a different distribution p(xf), the total entropy production is 

As = In p(xi) - In p(x f ) - /3Q (43) 

or, using equation (|38|) . 

As = \np(xi) - kip(xf) - U cq (x f ) + U cq (xi). (44) 

Thus, in the absence of an explicit time dependent driving force, the stochastic entropy 
is a state function which depends only on the initial and final states. Moreover, if we 
prepare the initial state in equilibrium p(xi) = P eq (xi) = Ne~ Ucci ^ Xi \ we immediately 
conclude that As = for each stochastic trajectory and for arbitrary value of a. The 
relevance of the potential U eq in the entropy production was recently pointed out in [38] 
for the Stratonovich prescription. Here, we generalize this concept to the general case 
of the a-prescription where the time reversed process is not trivial. It is interesting to 
note that a similar expression to ( f44l) . however quite cumbersome, could be obtained 
for a stationary out-of-equilibrium potential. The complication in this case is that the 
stationary potential depends on J, and the action does not transform as easily as in 
equation (1381) under time reversal, since J — > — J. 
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5. Correlation functions, Grassman variables and hidden symmetries 

In order to push forward our analysis of equilibrium properties it is convenient to have 
a formalism to compute correlations and response functions. However, the action of 
equation (l35j) and the difficulties introduced by the generalize Ito calculus (eq. (134]) ) 
make the functional formalism in x(t) quite cumbersome. For this reason it is convenient 
to introduce auxiliary stochastic variables (a function ip(t) and two Grassman functions 
£(£),£(£)) extending in this way the functional space of stochastic functions. The 
immediate consequence is a formalism that does not depend on the explicit value of 
a, simplifying its mathematical properties. As a by-product, hidden symmetries of the 
system can be put in evidence. These symmetries impose non-perturbative constraints 
on correlations functions and are related to equilibrium properties of the system like, 
for instance, the fluctuation-dissipation theorem. 

5.1. Generating functional 

We are interested in the calculation of correlations of functions of the stochastic variable 
at equal or different times. For instance, we could want to compute (A(x(t))), where 
A(x(t)) is an arbitrary function of x, or two point functions {A(x(t))B (x(t'))) . 
In the path integral approach these mean values can be written as 

(A(x(t))) = f dxidx f P(xi) f Vx det _1 (#) A(x(t)) e~ s[x] , (45) 

J Jx{U) = Xi 

where S[x] is given by equation (13T1) and U < t < tf. This equation takes a simpler 
form if we consider that the system is prepared in equilibrium at t — t^. In this case, 
P(xi) = P eg (xi). Assuming an ergodic evolution, we consider that 

Peq{Xi) = lim P(Xi,ti\X- T , -T) (46) 

for an arbitrary value of X-t- Therefore, using equation (l4"5l) and equation (|4"6"1) we can 
write the mean value as 

(A(x(t))) = (Vx det-\g) A(x(t)) e~ s[x] , (47) 

where we have extended the time integration from — oo to oo and, now, there are 
no more constraints in the functional integration. It is important to recall that this 
simpler expression is only possible because we sample the initial condition x(ti) with 
the equilibrium distribution P eq {xi). Provided this is satisfied, we can compute any 
mean value or correlation function from the generating functional 

Z(J) = Jvx &et-\g)e- S[x]+ ^ dt ' J{t ' )x{t,) , (48) 

where J(t) is a source with compact domain, that is, it adiabatically goes to zero away 
from an interval (ti,tf) where we will compute the correlation functions. 

In this context, any total derivative term in S[x] does not contribute to the dynamics 
of any observable. A total derivative term only contributes to a constant prefactor of a 
correlation function. In this way, a system described by equation ( HHl) is automatically 
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invariant under time reversal since S and S = TS differ just in a total time derivative 
term. On the other hand, if we impose the constraint x(+oo) = x{— oo), the action is 
truly time reversal invariant S[x] = S[x] = S. Again, it is important to observe that 
this happens because we have prepared the system at equilibrium at some initial time 
ti and we are computing correlation functions at later times t > U. For any other initial 
distribution we need to go back to equation (|4~5]) to compute mean values. 



5.2. Time reversal and entropy production in terms of Grassman variables 

Equation ( ]48i) can be used to study general correlation and response functions of a white- 
noise multiplicative stochastic process whose dynamics is driven by equation ([1]). The 
a-dependence of the final equilibrium distribution, the structure of S and the fact that 
for general a the usual rules of calculus do not apply, make this method cumbersome. 
Nevertheless, we have shown [22] that, introducing an auxiliary "bosonic" variable (p(t), 
and a couple of conjugate Grassman variables £(£),£(£), we can avoid this inconvenience 
by alternatively using the generating functional formalism in this extended space. The 
generating functional can be cast in the following form, 

Z[J] = J VxVxpViVl e -si*,v*&+JdtJ(tMt) ^ ( 49 ) 

where the "action" S is given by 

/+°° f _ H - 1 

-oo dt I ~' m Jt m + f'WWK® + 2 V{t)2g{x)2+ 
dx 



+ i<p(t) 



(jt -f(x) + g(x)g'(x)mm 



(50) 



This action in the extended functional space (x(t), ip(t), £(t), £(t)) is completely equi- 
valent to the usual formalism presented in §4.11 In fact, functional integrating over 
Grassman variables and over the variable tp we get the generating functional (I48p . with 
the action given by equation fl3T|) . The advantage of this extended formulation is that it 
does not depend explicitly of the a-parameter and then, the calculus rules are the usual 
ones for any a. Hence, all the complication of the stochastic calculus associated with 
the definition of the Wiener integral is now codified in the anti-commuting Grassman 
algebra, implying, for instance, that £(£) 2 = £(£) 2 = 0. Of course, the parameter a 
reappears when we need to properly define Grassman retarded Green's functions. In 
fact, equal-time correlations functions are ill-defined, forcing us to set (£(£)£(£)) = ol. 
One of the advantages of this procedure is that, while the Grassman variables are 
not integrated, we can work out any calculation without explicitly indicate an specific 
convention. This is particularly useful to study Markovian white-noise processes, where 
the equilibrium distribution depends on a. This situation is very different from non- 
Markov processes [28], in which (£(£)£(£)) = 0. 
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The time reversal transformation, equation ( 129ft . has a simpler form in terms of 
Grassman variables, 

' x{t) ->■ x(-t) 

T { m - v(-t) (51) 

where 77, 77 are the time reversed Grassman variables. There is a subtlety in these 
transformations which involves the extremal values of the Grassman variables. While 
lim^ioo £(t)£(t) = a, the time reversed variables satisfy lim^ioo fj(t)r)(t) = 1 — a. 
This fact has no importance at all when computing correlation functions. However, 
it is relevant to analyze entropy production and other relations that depend on total 
derivatives. 

The real advantage of the Grassman representation is that almost any important 
symmetry becomes a linear transformation. It is simple to check that, in this 
representation, T 2 = I and 

S - S = Uf q - U l eq + £^ In g 2 d - (&) dt. (52) 

Then, the entropy production As, associated with each trajectory, is given by 

Jng 2 -(tt) dt. (53) 

We see that, in the extended space of trajectories (x, (p, £, £), the entropy production 
associated with each trajectory is not zero. This is due to the apparently extra degrees 
of freedom introduced by the Grassman variables. On the other hand, upon functional 
integration we have that 

^<ft) = 0, (54, 

due to fermionic number conservation and the total entropy production is zero, as 
it should be in an equilibrium evolution. This opens the interesting possibility of 
interpreting an entropy production as a spontaneous breaking of fermionic phase 
(Gauge) invariance like, for instance, in superconducting phase transition. This effect 
cannot happen in our example of one stochastic variable, but in multiple variables 
systems [39], there are quartic interactions between Grassman variables opening this 
attractive possibility. 

5.3. Linear response 

It is interesting to compute the linear response of the system to an external perturbation. 
To see this, we slightly perturb the system out of equilibrium 

V{x) ->■ V{x) - h(t)x(t), (55) 

and compute the dynamic susceptibility 



xl,tj " 8h{t>) 



(56) 

h=0 
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Introducing (1551) into (1501) and computing equation (1551) we find 

X (t, t') = i(x(t),g(x')Mt')) ~ (x(t),g(x')g'(x>)at')m) (57) 
In the case of additive noise, where g = constant, the response function has the simpler 
form, x(tjt') ~ ( x (t)i V 9 (^'))- For this reason, the auxiliary function <p(t) is usually 
called the response variable. However, for multiplicative noise, the response is more 
complex since involves, not only correlations functions of g(x), but also contributions 
from the Grassman sector of the model. This is a direct consequence of multiplicative 
white-noise since a variation of the "external" potential is modifying the fluctuations 
properties of the system. This means that the susceptibility explicitly depends on a. 
From equation (l57l) . it is immediate to recognize a "natural response variable" 

m = g 2 (x) Lt) + i^fmm) , m 

in such a way that, 

X (t,lf)=i{x(t),i<f(t')) (59) 
We will show that (p(t) has a key role in the supersymmetric formulation of the model. 

5-4- Supersymmetry and covariant (superfield) representation 

The path integral formalism is useful to make evident symmetries of the stochastic 
process. For instance, the action given by eq. (1501) is invariant under the transformation, 

6x = M , = , (60) 

6£=i\(p , 6<p = 0, (61) 

where A is an anticommuting parameter. This nilpotent transformation (5 2 = 0) is the 
famous BRS [10] symmetry, discovered in the context of quantization of Gauge theories. 
In the present context, it simply enforces probability conservation, Z(0) = 1. 

There is another set of important symmetries related with equilibrium properties 
which is called supersymmetry. To display it explicitly, it is convenient to work with the 
natural response variable introduced in equation ( 1581) . Therefore, we make in equations 
( 149]) and (!50|) the funcional change of variables (x,(p,£,£) — > (x,<f,T],fj), defined by 

m = g 2 (x( t )) + i dln9 ^ {t)) mm) m 

v(t) = g(x(t))t(t) (63) 

m =g(x(t))l(t). (64) 

The last two transformations make the Jacobian trivial, T>(pT>^T>^ = TXpDrfDf], since 
T>Lp = det 2 (g)T>(p and the fermionic measure reads Vr]Vf] = det~ 2 (#)£>££>£. 
The generating functional is then given by 

Z[J] = J VxVtpVrfDf) e -s[x,v,v,v]+J dtJ(t)x(t) ^ ^ 
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where the action S[x, <p, r], rf\ reads 

r+°° ( d / \ 1 

s= I dt l-m^jtig v) + f'9~ 2 fiv + ^ 2 + 



DO 



+ ig V 



dx g _ 

-77 - f - —rm 

dt g' 



+ ^(i-f)w\, (66) 



where we have drop-off function's arguments just to simplify notation. 

At first sight, this action written in the space (x, <p, rj, ff) seems to be more complex 
than the preceding one. However, these new physical variables can be rearranged in 
order to explicitly display SUSY. 

First, it is convenient to write f(x) and g(x) in terms of two "potentials", V(x) 
and T(x), defined as follows 

/(*) = - (67) 

g-\x) = ^M. (68) 

Then, we collect the variables x, tp, fj, rj, in the definition of the superfield 

9, 9) = x(t) + 6r](t) + f}(t)9 + i(p{t)99. (69) 

where we have introduced two "temporal" Grassman variables 9 and 9. 

The generating functional, equation (166|) . can be re- written in terms of $(£, 9, 9) as 

Z[J] = J V® e~ 5[ * ]+ / dtd6d§ ; (70) 

where S[Q] depends on the potentials V($) and T($) by 

= J dtd9dd = J dtdddd IdT[$\DT[$\ + ^[$]| . (71) 

We have defined the covariant derivatives 
d Id d 

D = Te • D =2W- e m' < 72 » 

which satisfy D 2 = D 2 = and (d, Dj = — ^. The action (JTIjl . when written in 
components form, is completely equivalent to equation f loU]) or, after integration of the 
Grassman variables, to equation ( 136|) . as it is shown in the Appendix. 

It is now immediate to verify that (j7ip is invariant under transformations of the 
supersymmetry group [26], whose generators are 

Q = d m • «=y* +s i • {«■«}=!■ < 73 » 

The symplicity of the superfield representation resides in the fact that each SUSY 
transformation is a translation of the temporal variables (t, 9, 9). 

Then, the invariance of the last term of equation (1711) is trivial, since U($) does 
not depend explicitly on (t,9,9). To demonstrate the invariance of the kinetic part of 
the action, it is useful to have in mind the graded algebra, 

{Q, D} = {Q, D) = {Q, D) = {Q, D) = 0, (74) 

Q 2 = Q 2 = . (75) 
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Thus, 

5 {DT[$]DT[$]} = D5T[$]DT[$] + DT[$]DST[$]. (76) 

On the other hand, 5T[<&] = eQpr[$>], where e is an arbitrary Grassman parameter and 
Qp may represent Q or Q defined in equation (173]) . In this way, 

5 {DT[$]DT[$]} = e {-D {Q P T[<$>}) DT[<$>] + DT[$]D (QpT[$])} (77) 

Finally, using the algebra ( 1T4"1) . we find, 

5{DT[$]DT[$]} = eQp{DT[$]DT[$]}, (78) 

which is a total derivative term, demonstrating in this way, the invariance of S. 

For additive processes, the diffusive potential is linear, T(x) ~ x, and eq.( ITT|) 
reduces to the usual action defined with a single superpotential V(<&). In this case, 
the tadpole theorem [22] guarantees that the stochastic evolution does not depend on 
a. However, multiplicative processes (non-linear T) induce derivative couplings in the 
superfield. These couplings are responsible for the a-dependent evolution which leads 
to the equilibrium distribution of equation (fl5|) . 

5.5. SUSY Ward-Takahashi identities and equilibrium 

The study of the symmetries of the action enables us to deduce properties of 
equilibrium dynamics of the system. Each of the three generators of SUSY, {Q, Q, d t }, 
imposes several non-perturbative constraints on correlation functions. To obtain these 
constraints, we will use the Ward-Takahashi identities 

($(*!, 0i, 0i)$(* 2) e 2 , e 2 )) = (e eG e u h)e eG $(t 2 , 02, e 2 )), (79) 

where e is an arbitrary parameter and G is one of the generators of SUSY. Note that, 
if G is Q or Q, then e is a Grassman parameter, while, if G = d t , e is a commutative 
parameter. 

The Ward-Takahashi identity related with the generator d t reads, 

($(*!, e u ^)$(t 2 , e 2 , e 2 )) = Wi + e,e 1 , + e, e 2 , e 2 )), (so) 

implying time-translation invariance as 

+ d t2 )($(t u e h hMh, e 2 , e 2 )) = 0. (si) 

Then, any two-point correlation function depends on time differences t 2 — t\. This is a 
necessary condition for an equilibrium evolution. 

On the other hand, the generator of 9 translations, Q, induces the identity 

($(*!, 6 1 , 0i)$(* 2 , 6 2 , h)) = Wi, Oi, 0i + e)<S>(h + e, 6 2 , 2 + e)>, (82) 

that, when written in components, imposes the relations 

(x(ti)itp(h)) = (f)(t 2 )v(ti)), (83) 

mmh)) =o. (84) 
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All components, involving only one of the Grassman variables, r] and f), vanish because of 
fermionic number conservation. The first equation relates the physical linear response 
of the system with the Grassman two-point correlation function. Since the physical 
response is causal, i. e. , {x{ti)i(p{t 2 )) = for t\ < t 2 , (see equation (|59|) ). we are 
forced to choose the retarded prescription in computing Green's functions of Grassman 
variables. On the other hand, the equal-time Green function is not well defined, forcing 
us to adopt a prescription (fj(t)r)(t)) = ag 2 ^^), where arbitrary initial value. 

This presciption, when written in the original variables, is (£,{t)l;(t)) = a, allowing us 
to identify a with the needed prescription to define the Wiener integral. Notice that, in 
the case of non-Markovian processes, (£(£)£(£)) = 0, without any ambiguity. 

Finally, the invariance of the action generated by Q, which arise from the identity 

($(tl,01,0l)$(*2,M 2 )) = 

+ 0x + e/2, hmt 2 + e6 2 , 6 2 + e/2, 2 )>, (85) 
results in the non-perturbative constraints, 

(d tl x(ti)x{t 2 )) = \(vp{ti)x(t2)) + l(v{ti)fj(t 2 )), (86) 
(d tl x(timh)) = (V(ti)d t2 fj{t2)). (87) 

This constraints imply another important equilibrium property, the fluctuation-dissipa- 
tion theorem, which relates the spontaneous fluctuations of the system with its response 
to an external perturbation. To explicitly see this, we substitute relation (l8~3l) into ( 1861) . 
obtaining 

d t (x(t)x(t')} = \mt)x{t')) - \{x{tmt'))- (88) 

Using the equilibrium time translations invariance of the correlation function C(t,t') = 
(x(t)x(t')), this identity becomes 

\{d t - d t ,)C(t - = ~ X (i' - t) - X (t - t'), (89) 
which, bearing causality in mind, can be rewritten as 

x(t - f) = -{d t - d t ,)c(t - f )e(* - f), (90) 

where 0(t) is the Heaviside step function. Equation fl90l) is a traditional form of the 
fluctuation-dissipation theorem. A more commonly used form of equation ( 19"0"|) is given 
in terms of the Fourier transform of the response function, 

(x(u)x(-u)) = (91) 

CO 

Then, the fluctuation-disipation theorem can be understood as a direct consequence of 
supersymmetry. We have shown that, even thought the correlation functions and the 
responses depend on a in a multiplicative process, the fluctuation-dissipation theorem is 
satisfied for any value of a, provided we correctly identify the equilibrium distribution, 
equation (fT5"]) and the time reversal transformation, equation ff29|) . These physical 
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concepts are mathematically codified in the invariance under supersymmetry trans- 
formations, that we have explicitly shown in this section. 

It is very interesting to note that extensions of the fluctuation-dissipation theorem, 
for out-of-equilibrium steady states, can be formulated A natural arising question 
is whether a related SUSY exists for these out-of-equilibrium states. 

6. Discussion and conclusions 

We have presented a detailed study of equilibrium properties of a Markovian multiplica- 
tive white-noise stochastic process. The stochastic process was modeled in its simplest 
version of a single stochastic variable, by means of the Langevin equation ([1]), defined 
by a drift term f(x) and a diffusion function g(x). We have completed the definition of 
the stochastic differential equation by using the Generalized Stratonovich Convention, 
characterized by a parameter < a < 1. 

The asymptotic equilibrium distribution was found by solving the stationary 
Fokker-Planck equation, imposing an asymptotic zero current probability. This na- 
turally conduce to an equilibrium potential U eq (x) (Eq. ffT5]) ) which explicitly depends 
on a. For a = 1 (Hanggi-Klimontovich interpretation), the potential coincides with 
the classical deterministic one, V(x), in agreement with the Boltzmann distribution 
for thermodynamical equilibrium. However, for any other convention, the equilibrium 
potential is more general, including, for instance, power-law distributions. Thus, the 
definition of equilibrium in the stochastic dynamical sense not necessarily coincides 
with the thermodynamical equilibrium concept. We have shown that, even for this 
generalized definition of equilibrium, the system satisfies the usual equilibrium properties 
such as detailed balance, microscopic reversibility and the fluctuation-dissipation 
theorem. 

Time reversal properties of the stochastic process are quite interesting. We have 
shown that, if the forward stochastic trajectory evolves with a definite value of a, the 
time reversed trajectory evolves with the conjugated prescription (1 — a). Therefore, in 
order to have a unique equilibrium distribution, the definition of time reversed stochastic 
process is given by the transformation fl29|) . where we need to change not only the 
velocity sign, but also the prescription and the drift force. 

We have shown that, using the equilibrium potential U eq and the time reversal 
transformation of equation ( |29|) . the stochastic process satisfies detailed balance relations 
for any value of the parameter a. It is convenient to note that, in equation ( 140L 
the forward transition probability is not equal to the backward one, except for the 
Stratonovich convention a = 1/2, where forward and backward trajectories coincide. 
We have also shown that, if the initial and final states of a finite time evolution are 
prepared in equilibrium, the entropy production of each stochastic trajectory vanishes 
for arbitrary prescriptions, verifying the related principle of microscopic reversibility. 
In order to show these properties, we have used the path integral representation of 
the stochastic process, where the central object is the hole stochastic trajectory. Also, 
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we were forced to intensively use an a-generalized Ito calculus, summarized in the en- 
dependent "chain rule" of equation fl34|) . 

Related with equilibrium properties, there are hidden symmetries in the stochastic 
process that can be better analyzed in the path integral approach, formulated in 
an extended functional space composed of commutative as well as anti-commutative 
Grassman variables. The advantage of this formulation is that it is prescription 
independent. That is, ambiguities in the definition of the multiplicative stochastic 
processes can be understood, not as a continuum limit ambiguity, but as a necessity to 
define equal time Grassman Green's functions. Then, provided the Grassman variables 
are not integrated out, it is not necessary to fix any prescription. In this way, we can 
use ordinary calculus rules, but instead we need to work with anti-commuting variables. 
In some sense, all the complexity of the Ito calculus is codified in the Grassman algebra 
in this formulation. 

In this context, we have computed the entropy for each stochastic trajectory (in 
the extended space) and we have shown that it is not zero, but proportional to a time 
derivative of the Grassman variables density (equation fl53|) ). This is due to the extra 
degrees of freedom we have added to extend the functional space. However, integrating 
over the Grassman variables, we immediately recognize that the entropy production 
is zero, as it should be for an equilibrium evolution. Interestingly, this cancellation is 
due to "fermionic number conservation" , in other words, due to the invariance under 
phase transformations. This allows us to speculate about the very attractive possibility 
of interpreting an out-of-equilibrium evolution, i. e. , a non-zero entropy production, 
as a consequence of a gauge symmetry breaking, very similar with a superconducting 
transition. Of course, this cannot happen in our present example, since we only have 
one stochastic variable and, consequently, the Grassman variables are not-interacting. 
Nevertheless, in multiple variable systems, quartic terms in the Grassman variables are 
certainly possible. 

Another interesting symmetry which appears in this formulation is supersymmetry 
(SUSY), expressed as a linear transformation, mixing commuting as well as anti- 
commuting variables. We have presented a supersymmetric covariant representation 
of the stochastic process. While this representation was known for additive processes 
and multiplicative non-Markov processes, the present case of Markovian multiplicative 
white noise was quite difficult to treat due to the different prescriptions available to 
define it, connected by a time reversal transformation. In this representation, the non- 
trivial time reversal transformation is represented by a simple linear transformation. 
The difficulty of the time conjugated prescriptions are confined to boundary conditions, 
irrelevant in the computation of correlation functions. 

Using this SUSY covariant formalism we have analyzed two-point Ward-Takahashi 
identities. One of them relates linear responses and fluctuations. We have shown that, 
while the linear response function and fluctuations are a-dependent, they satisfy the 
fluctuation-dissipation theorem for arbitrary values of the prescription a. Therefore, 
similar to the additive case, supersymmetry is a consequence of equilibrium evolution, 
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even thought we have different equilibrium evolutions for different prescriptions. 

Summarizing, we have presented a compact path integral formalism to deal 
with Markovian multiplicative white-noise systems, independently of the prescrip- 
tion used to define the Wiener integral. In its supersymmetric covariant form, 
it automatically encodes equilibrium properties. Therefore, it is very appropriate 
to organize perturbative calculations, since supersymmetry imposes non-perturbative 
constraints that should be verified at any order of perturbation theory. On the other 
hand, keeping in checked equilibrium properties, we can safely go forward in the study 
of out-of-equilibrium fluctuation relations described by this type of stochastic processes. 
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Appendix A. 

In this appendix we explicitly show that the action, written in the super-field formalism, 
equation flTTl) . is in fact equivalent to the time reversal invariant action S, given by 
equation ( f36l) . 

Let us define the scalar super-field 

$(t, e, 9) = x(t) + en(t) + fj(t)e + i(p{t)ee. (A.i) 

and the super-potentials V($) and T($), that we expand in powers of the temporal 
Grassman variables 9 and 9 as, 

r($) = r(a;) + r'(a;) + fj(t)e} + do {i<p(t)V(x) - v'\x)^} , (A.2) 

and 

V(<£>) = V(x) + V\x) [9r](t) + fj(t)0} + 99 {i(p{t)V'{x) - V"(x)f}r]} .(A.3) 
The covariant derivatives (1721) acting on the super-potential T($) take the form 

DT($) = -T'(x)fj - 9 {i0(t)T'(x) - T"(x)fjr]} (A.4) 

and 

+ 99 j t (V'(x) V ) . (A.5) 



The system Lagrangian is given by the 99 components in the form 

1 
2 



DTDT+-V} . (A.6) 
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Then, using equations (IA.3j) . (IA.4|) and fl A. 5[) . 

Dtdt 



- fjr'^- (r'ri) + up Jr'^-r + r'n^l - 



dt { dt 

- r'jv m + iy5 2 r 2 (a.?) 

and 

V^)\ Be ^i0(t)V(x)-V"(x)fjrj. (A.8) 
Replacing these expressions into equation (1A.6I) we find 

i --^( r, ') + ^H r+ ^' +r ' r "^}- 

■ r "it rm + ¥ 2r ' 2 ~ (A - 9) 

We make the following change of variables in the functional integral: first, we shift the 
variable <p, 

(pit) ^ Vl {t) - i2^- (A.10) 

whose Jacobian is one, since it is just a translation. Then, we rescale (ipi,rj,f)) in the 
following way, 

r)(t) -». (T'ix))- 1 ^) (A.ll) 

m -^mn^r 1 (A.12) 

^(i)^^)^))- 2 . (A.13) 

The Jacobian of the first two transformations is det 2 (T') and the Jacobian of the 
last transformation is det _2 (r'), then the two Jacobians cancel each other and the 
Lagrangian in the new variables reads, 

L — + up I — -r + k } - 

^ 2T' 2 *\Fdt 2T' 2 r /3 ^J 

Considering that, by definition, T'(x) = l/g(x) and, thus, T"{x) = —g'(x)/g 2 (x), 
we obtain 

L = - eC + \gV +i<f{x + \g 2 V + gg'tf} - 

l/ 9 ,\/- ldlnq 2 - Id , . . 

The last term is a total derivative. This term appears making a trivial transformation 
like t — y t + 1/299 in the V($) term. This is because a supersymetric Lagrangian is 
invariant up to a total derivative under SUSY transformations. This is a general feature 
of SUSY. A transformed supersymmetric action is invariat up to boundary terms. 
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Rewriting this last expression in terms of the original functions / and g we finally 

get 

1 1 d 

L = 2#V+^ {i~f + gg'ti} + ttf'-2j t { V + &^9 2 } -(A.16) 

Upon integration and realizing that the last term is a total derivative one, we find 
Ssusy = dtL = S- l -(V + a \ng 2 )^ = S. (A.17) 

Where S is the action of equation (150]) and, then, S is equal to equation ( 136|) . as can be 
seen comparing (1A.17j) with (155]) . as we want to demonstrate. 

To compute correlation functions the range of integration is infinite and S = S. 
However, for finite time we see that the supersymmetric action is equivalent to the time 
reversal invariant action S. 
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